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While NMR measurements of nuclear energy spectra are routinely used to characterize the static 
properties of quantum magnets, the dynamical information locked in NMR 1/Ti relaxation rates 
remains notoriously difficult to interpret. The difficulty arises from the fact that information about 
all possible low-energy spin excitations of the electrons, and their coupling to the nuclear moments, is 
folded into a single number, 1/Ti. Here we develop a quantitative theory of the NMR 1/Ti relaxation 
rate in a collinear antiferromagnet, focusing on the specific example of BaFe2As2. One of the most 
striking features of magnetism in BaFe2As2 is a strong dependence of 1/Ti on the orientation of 
the applied magnetic field. By careful analysis of the coupling between the nuclear and electronic 
moments, we show how this anisotropy arises from the "filtering" of spin fluctuations by the form- 
factor for transferred hyperfine interactions. This allows us to make convincing, quantitative, fits 
to experimental 1/Ti data for BaFe2As2, for different field orientations. We go on to show how 
a quantitative, angle-dependent theory for the relaxation rate leads to new ways of measuring the 
dynamical parameters of magnetic systems, in particular the spin wave velocities. 

PACS numbers: 67.80.dk 76.60.-k 76.60.Es 



I. INTRODUCTION 

Nuclear magnetic resonance (NMR) has a long history 
as an experimental probe, with numerous uses through- 
out physics, chemistry and medicine. More than 60 years 
after its discover^i^, it remains one of the most powerful 
techniques for investigating solid state systems. NMR 
spectra measurements of the nuclear energy level split- 
ting are well understood, and provide a wealth of quan- 
titative information concerning the static properties of 
magnetic materials^. However, the dynamic properties 
of these materials are more difficult to access. While 
much of the relevant information is encoded in the NMR 
relaxation rate 1/Ti, it can prove difficult to extract. 

The problem with interpreting measurements of 1/Ti 
is that all possible fluctuations of the electron moments, 
as well as the details of the coupling to the nuclear mo- 
ment, are folded into a single number. Despite many 
decades of stud}*^— , and some notable successes, the- 
ory has thus far not developed the level of sophistication 
required to fully utilise the information stored in these 
measurements. In this paper we address this problem by 
developing a quantitative theory of the 1/Ti relaxation 
rate in the magnetically ordered phase of BaFe2As2. 

The Fe-pnictide materials in general, and BaFe2As2 in 
particular, have been a field of much activity in recent 
yearsi^Ti^U. In BaFe2As2 layers of FeAs alternate with 
planes of Ba. While at room temperature the undoped 
materials support a tetragonal /4/mmm, paramagnetic 
phase, at approximately 135K there is a structural dis- 
tortion to an orthorhombic Fmmm phase, and subse- 
quently a magnetic ordering to form a striped, collinear, 
antiferromagnelii^. Chemical doping, such as the sub- 
stitution of Co for Fe, suppresses the magnetism of the 
parent compounds and leads to an intriguing supercon- 
ducting state^l. 

One of the many puzzling features of magnetism in 
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FIG. 1: (Color online). Experimental measurements of the ^^As 
1/Ti relaxation rate in BaFe2As2, as reported in Ref. [l2|. Relax- 
ation rates were measured with external magnetic field applied both 
perpendicular (squares) and parallel (circles) to the FeAs planes. 
Except at the lowest temperatures [inset] , the relaxation rate mea- 
sured for field parallel to [001] is significantly lower than that for 
fields parallel to [110], and has a qualitatively different temperature 
dependence. 



BaFe2As2 is a significant reduction in the ^^As 1/Ti 
NIVIR relaxation rate when the external magnetic field is 
applied perpendicular, rather than parallel, to the FeAs 
planes [Kitagawa et ali^^ reproduced in Fig. ([1])]. This 
sensitivity of 1/Ti to the orientation of magnetic field 
cannot be explained by any existing theory of NIVIR re- 
laxation rates. 

In this paper we develop a quantitative theory of the 
^^As NIVIR 1/Ti relaxation rate in the magnetically or- 
dered phase of BaFe2As2, building on the earlier ideas of 
IVtoriya^^, and IVtila and Rice^*^. We find that the tensor 
nature of transferred hyperfine interactions between elec- 
tronic and nuclear spins leads to a "filtering" of the spin 
fluctuations seen at the ^^As site, which in turn depends 
on the orientation of the magnetic field used in NIVIR 
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experiments. Consequently, the NMR 1/Ti relaxation 
rate has a qualitatively different temperature dependence 
for fields applied parallel and perpendicular to the FeAs 
planes. This theory is developed in absolute units, and 
its predictions are found to be in quantitative agreement 
with experiment, for both orientations of the magnetic 
field. 

In this context, it becomes meaningful to talk about 
"angle-resolved" measurements of 1/Ti, and we there- 
fore generalize our theory of ^^As NMR 1/Ti to treat ar- 
bitrary orientation and magnitude of external magnetic 
field. We use this more general theory to make specific, 
quantitative, predictions for the shape of the 1/Ti surface 
as a function of field orientation at fixed temperature, 
and the field-dependence of the 1/Ti relaxation rate as a 
function of field strength, at fixed field orientation. 

The combination of angular resolution and absolute 
units also permits quantitative information about dy- 
namical properties of the system to be extracted directly 
from 1/Ti measurements. As an illustration, we show 
how angle-resolved ^^As NMR 1/Ti measurements could 
be used to measure all three components of the spin wave 
velocity in BaFe2As2. 

While this paper is concerned with the specific example 
of ^^As NMR in BaFe2As2, many of these results general- 
ize straightforwardly to other collinear antiferromagnets, 
and corresponding results can be derived for more exotic 
magnetic states. As such, angle-resolved measurement of 
NMR relaxation rates promises to be a powerful probe 
of both conventional and unconventional magnetism. 

The paper is structured as follows : In Section |TT1 we 
present some of the basic facts about antiferromagnetism 
in BaFe2As2, and briefly review existing theories of NMR 
1/Ti rates in antiferromagnets. In Section HITl we intro- 
duce the idea of angular resolution in 1/Ti measurement, 
and develop a theory for the ^^As 1/Ti relaxation rate 
in the magnetically ordered phase of BaFe2As2, for fields 
applied in both [110] and [001] directions. In Section HVl 
we show how this theory can be extended to treat ar- 
bitrary field orientation and strength. In Section |Vl we 
propose a scheme for measuring spin wave velocities di- 
rectly from NMR 1/Ti relaxation rates, based on the 
results of Section IIVI We conclude in Section IVII with 
a discussion of some of the wider implications of these 
results. 



II. BACKGROUND TO NMR 
MEASUREMENTS ON BAFE2AS2 

A. Low temperature magnetism in BaFe2As2 

Before delving into the details of the NMR re- 
laxation rate, we briefly review the nature of the 
low temperature magnetic state in BaFe2As2. Neu- 
tron scattering measurement^— reveal a commen- 
surate, collinear antiferromagnetic ground state be- 
low 135K, with ordered moment mo ~ 0.87/iB [3 5 



and ordering vector Q = (yr/ao, 0, tt/cq)^. [Here 
(ao, bo, Co) = (2.80 A, 2.79 A, 6.47 A) are the lattice con- 
stants of the orthorhombic lattice of Fe atoms within 
BaFe2As2]. The ordered magnetic moment lies along the 
crystallographic a-axis, as shown in Fig. (j2j). 

A single branch of low-energy spin wave excitations 
with dispersion, 

fi^q = /a2+ J2 ^lida-Qcr, (1) 

Y a=a,b,c 

is found above an anisotropy gap A ^ 9.SmeV [l6|. The 
spin wave velocities are anisotropic with Va ^ Vb Vc 
These spin-wave excitations become diffuse at 
higher energies, and merge into a broad continuum of 
incoherent spin excitation^. 




FIG. 2: (Color online). Below 135K BaFe2As2 exhibits 
collinear antiferromagnet order with characteristic wave vector 
Q = (7r/ao, 0, tt/co). Within this ordered state Fe moments — in- 
dicated here by red, horizontal arrows — are orientated along the 
crystalographic a-axis. The effective magnetic field induced at the 
As nucleus by these Fe moments, hint (blue, vertical arrows), is ori- 
entated along the crystalographic c-axis, and alternates in direction 
between the As sites. 

Constructing a theory of NMR relaxation rates is com- 
plicated by the fact that BaFe2As2 is both an antifer- 
romagnet with a sizeable ordered moment, and a metal. 
The existence of a Fermi surface in BaFe2 As2 implies that 
it must support gapless particle- hole pairs, as well as the 
coherent spin-wave excitations seen in neutron scatter- 
ing. These incoherent particle-hole pairs will contribute 
to the 1/Ti relaxation rate in their own right and, a pri- 
ori, might be expected to couple strongly to spin waves. 

In Section IIIII we model the magnetic excitations of 
BaFe2As2 using a low temperature field theory written 
in terms of the hydrodynamic parameters Va^ v^, Vc^ A 
and the transverse susceptibility x±- This field theory 
respects all of the symmetries of the magnetic ground 
state, correctly reproduces the low temperature disper- 
sion, Eq. ([1]), and provides self-consistent predictions for 
all related magnetic properties. As such it offers a cor- 
rect low-energy effective theory of spin-wave excitations 
in BaFe2As2, regardless of the microscopic details of the 
material. 

The effect of incoherent particle- hole pairs, neglected 
in this theory, is most evident in the linear-T behavior 
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of 1/Ti at temperatures T < 50K [cf. inset to Fig. [T]. 
We have previously argued that this contribution to 1/Ti 
remains hnear in T at higher temperatures, and can be 
safely fitted at low temperatures and subtracted from the 
datapi^. We return to this argument below. 



B. Introduction to NMR relaxation rates in 
antiferromagnets 

NMR measurements of the average internal magnetic 
field are well understood, and angular resolution is rou- 
tinely used to determine the static properties of magnetic 
materials^. The technique relies on the Zeeman splitting 
of the nuclear energy levels by the effective magnetic field 
at the nuclear site, hpuc, as described by. 



= 7ArM.(hnuc) 



(2) 



Here I is the nuclear moment, /+(~) is a nuclear spin 
raising (lowering) operator, and 7Ar is the gyromagnetic 
ratio of the nucleus in question. 

The effective magnetic field experienced by the nuclear 
moment hnuc = hgxt + hint, is the vector sum of the ex- 
ternal magnetic field applied to the sample, hgxt, and an 
effective internal magnetic field, hint- This internal field 
encapsulates the effect of interactions between the nu- 
clear moment and the surrounding electrons. However, 
while the external field hgxt can be considered to be static, 
hint fiuctuates on a time scale set by the electrons, and 
so contributes to relaxation of the nuclear spins. 

In general, three different types of interaction can con- 
tribute to the effective internal field hint- Firstly, for 
the nuclei of magnetic atoms, there is an on-site hyper- 
fine interaction. Secondly, for the nuclei of non-magnetic 
atoms, there is a transferred hyperfine coupling between 
the nuclear moment and the spin of neighbouring elec- 
trons. Finally, we can also consider the dipolar interac- 
tion between nuclear and electronic spins. This is weak 
for small values of electronic and/or nuclear spin, but 
long-ranged. In the case of NMR in BaFe2As2, we focus 
in particular on the transferred hyperfine coupling be- 
tween an ^^As nucleus, and the electrons of the four Fe 
atoms which surround it. 

In NMR 1/Ti measurements, the population of the 
Zeeman split nuclear energy levels is driven out of equi- 
librium by a radio- frequency pulse. These nuclear spins 
then return to equilibrium over a characteristic timescale 
Ti, which is set by their interaction with electrons — 
specifically by the transverse fields /i+ and h~ in Eq. (j2j). 
We characterize this process by the standard 1/Ti = 2W 
relaxation rate, as defined in |22l-[2^. 

The quest for a theory of NMR 1/Ti relaxation rates 
in antiferromagnetic materials now spans almost six 
decades. Moriya, writing in 1956, was the first to re- 
alise how the Raman scattering of antiferromagnetic spin 



waves from nuclear moments can lead to a finite nuclear 
spin relaxation rate^— . Beeman and Pincus later ex- 
tended Moriya's semi-classical theory to include quan- 
tum fiuctuation effects^. The next major breakthrough 
was due to Mila and Rice, who realised that the in- 
direct coupling between electronic spins and the nuclei 
of non-magnetic atoms can act as a "filter" on spin 
fiuctuations^'^. This made it possible to understand for 
the first time why the 1/Ti relaxation rates of different 
nuclei in the same compound can have qualitatively dif- 
ferent leading temperature dependences. 

The theory of NMR relaxation rate presented in this 
paper extends Mila and Rice's idea of "filtering" of spin 
fiuctuations by showing how the action of the "filter" is 
strongly dependent on the orientation of the externally 
applied magnetic field. We develop this theory in abso- 
lute units, which allows us to make quantitative predic- 
tions for comparison with experiment. Using this theory, 
we are able explain the temperature dependence of ^^As 
relaxation rates in BaFe2As2, for fields applied in both 
the [110] and [001] directions'^. This task is simplified 
by the very large anisotropy gap A lOOi^T in the spin 
wave spectrum of BaFe2As2, which forbids additional re- 
laxation processes arising from the absorption or emission 
of spin waves^ii. 

We note that a theory of NMR 1/Ti relaxation rates in 
magnetic Fe pnictides has also been advanced by Ong et 
al^. However the theory we present goes much further 
and allows quantitative comparison with experiment^. 

We also remark that BaFe2As2 is not the first material 
in which the NMR 1 /Ti relaxation rate has been found to 
depend on the direction in which the magnetic field was 
applied — although it is, to the best of our knowledge, 
the first antiferromagnet. Angle dependence of 1/Ti re- 
laxation rates has previously been observed in ^^Cu NMR 
of the cuprates YBa2Cu307 and YBa2Cu408, within 
their low-temperature superconducting stat o^^i^^ . The- 
len. Pines and LvM have suggested that the anisotropy 
in 1/Ti follows from an anisotropy in the on-site cou- 
pling between the electron and nuclear moments of ^^Cu. 
In this paper we focus on a non-magnetic ion, ^^As, for 
which there is no on-site coupling, and explain the an- 
gle dependence in terms of transferred hyperfine interac- 
tions. 



III. QUANTITATIVE THEORY OF ^^AS 1/Ti 
FOR FIELD PARALLEL TO [110] AND [001] 

We now develop a theory of 1/Ti relaxation rates in 
BaFe2As2, with the specific goal of obtaining quantitative 
fits to high-quality ^^As NMR data for magnetic fields 
parallel to both the [110] and [001] directions'^. We fo- 
cus on results for the relaxation rate at low temperatures 
and in the magnetically ordered phase, and follow a rela- 
tively direct path to the results needed to compare with 
experiment. The development of a more general theory 
for arbitrary field orientation is postponed to Section HVl 
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If the nuclear field, hp 



is orientated in the z- 



direction, then time dependent perturbation theory leads 



to^, 



1 



Ti(hext) 



+ (WntW.C(0)})], 



(3) 



where huoo is the splitting of the nuclear energy levels and 
hint is the effective field at the nuclear site due to inter- 
action with the surrounding electron moment^. Only 
the X and y components enter the expression, since they 
are the ones that couple to the nuclear spin raising and 
lowering operators, in Eq. (|2j). 

The internal field, hjnt, can be re-expressed in terms of 
the electronic degrees of freedom as. 



hint(t) = ^A.mi(t), 



(4) 



where i sums over the electron moments (m^) that couple 
to the nuclear spin, and A. \T / (ib] is a rank two tensor 
describing the coupling between the nuclear and electron 
moments. We refer to ^.as the nuclear-electron coupling 
tensor, and note that its components. 



A. = 



ill Al2 Al3 
A21 A22 A23 

"^^i 

^31 ^32 ^33 



can be measured by Knight shift experiment^. The full 
interaction Hamiltonian for a nuclear moment, I, is given 

by, 



(6) 



Starting from Eq. (|3]), we now show how to derive 
a working expression for the relaxation rate by making 
use of the fluctuation dissipation theorem. Eq. ^ and 
Eq. (j3j) can be combined to give. 



7^ 



Ti(hext) 2 



/ 



dt 



N 



{(A.mq(t))-,(A.m_q)-} 
{(A.mq(^))^(A.m_q)n 



(7) 



where the electron magnetic moments have been Fourier 
transformed using, 



mq(t). 



(8) 



Fluctuations of the electronic moments can be charac- 
terised by introducing a dynamical structure factor. 



5«'^(q,u;o) = ^dt e-«* ({m« (t), m^j) . 



(9) 




FIG. 3: (Color online). The local environment of the As atom in 
BaFe2As2. The As atom (purple) experiences an average local field 
(hnuc) (blue arrow) that arises from a combination of an externally 
applied field hgxt, shown here in the a-direction, and an internal 
field (hint). The internal field points on average in the c-direction, 
and is due to the interaction with the electron magnetic moments 
(average position shown by red arrows). After being excited by 
a radiation pulse, the As nuclear moment I (green arrow) relaxes 
back towards alignment with the magnetic field. 



(5) It follows that the relaxation rate is given by. 



(10) 

where ^, = {x, z} and we have defined, 

A^^ = Y,e^-^^Af. (11) 

The fluctuation-dissipation theorem relates the struc- 
ture factor of Eq. ([TQ|) to the dynamic susceptibility, de- 
fined by. 



X 



(12) 



This allows the relaxation rate to be rewritten as, 
1 



^{x^^(q,c^o)} 



(13) 



where, in taking the limit cjq ^ 0, we have assumed 
that the energy of the nuclear transitions is negligible 
compared to the typical spin wave energies. 

Since both the susceptibility, Eq. (fT2]) . and the interac- 
tion between the nuclear and electron moments, Eq. (j5j), 
are tensors, it is important to keep track of the coordi- 
nate basis in which they are represented. In the above 
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expressions, both are represented in the {x^y^z) coor- 
dinate system, which is defined by ahgning the 2:-axis 
with hnuc, the nuclear magnetic field. However, they are 
most simply measured and calculated in the (a, 6, c) co- 
ordinate system, which is aligned with the crystal axes. 
Fig. (j3j) shows the orientation of these two coordinate 
systems when hgxt is parallel to the a-axis. In Section HVl 
we consider an arbitrary direction of hgxt, and introduce 
rotation matrices in order to map between the {x^y^z) 
and (a, 6, c) bases. Here we specialise to the two cases 
measured by Kitagawa et ah hgxt orientated in the [110] 
direction; and hext parallel to the c-axis. 

In what follows, we make the assumption 
I hext I ^ I hint I . This means that hpuc ^ hgxt and 
the z-axis is aligned with the external magnetic field. 
The discussion of an arbitrary magnitude of external 
field is postponed to Section |IVl If hext is applied 
in the c-direction then the (x, z) and (a, 6, c) coor- 
dinate systems are equivalent. In consequence, the 
tensors appearing in Eq. ([13]) can be transformed 
into the (a, 6, c) coordinate system by the substitution 
{x ^ a^y ^ z ^ c). If instead hext is applied along 
the [110]-direction, a set of rotation matrices is required 
to relate the tensors expressed in the two coordinate 
systems. Again we postpone the details of this to 
Section IIVI but, at a schematic level, we make the 
substitution {x ^ {h — a)/ >/2, y ^ c^z ^ (a + 6) / ^/2) . 

We have already argued that, due to the sizeable gap in 
the spin-wave spectrum, the relaxation is dominated by 
the scattering of spin wave excitations. This corresponds 
to picking out the longitudinal component of the suscep- 
tibility tensor in the crystallographic (a, 6, c) coordinate 
system. For collinear magnetic order in the a-direction 
this is the component X^^(q?^o) = XiKq^^o)- Thus the 
relaxation rate can be expressed as. 



moments. These two quantities can be developed inde- 
pendently, and then recombined to find the 1/Ti relax- 
ation rate. 



A. Form factor 

In this subsection we determine the form factors nec- 
essary to explain the experimental data shown in Fig. ([1]) 
[l2[. This requires an appreciation of the symmetry of 
the nuclear environment, and below we outline an analy- 
sis of the nuclear-electron coupling tensor similar to that 
carried out in [l2| . 

For the case of the As atom in BaFe2As2, the dominant 
interactioni^ with the electron system is via a transferred 
hyperfine coupling with the four nearest neighbour Fe 
electron moments, shown in Fig. (j3j). 

The environment of the As nuclear moment is thus 
invariant under the symmetry operations of the point 
group C2v' Referring to the labelling of the Fe moments 
shown in Fig. (|3]), the nuclear-electron coupling tensor 
for the first Fe site can be written as. 



j^aa J^CLb ^cic 
J^ba J[bb J^c 
J\ca j^h j^c 



(17) 



Refiection symmetry in the 6c-plane allows the tensor for 
the second Fe site to be determined as, 



4 



J^Ob J^c 

J^Ci j^c 



(18) 



By refiection in the ac-plane. 



ujo^o2N ^ 

qGPMBZ 

^^{xii(q, ^o)} 



Ti(hext) 



J'(q, hext)- 



(14) 



where J-'(q, hext) is a form factor for the nuclear-electron 
interaction, which acts as the "filter" of spin fiuctuations. 
The sum is over all q vectors in the Fe-paramagnetic 
Brillouin Zone (PMBZ). For external field in the [110]- 
direction. 



\ba\ ( Aaa _ 



Aca Aca 



while for external field in the c-direction, 

TT/p, hc \ _ Aaa Aaa , Aba Aba 



(15) 



(16) 



Eq. (p!4|) gives the relaxation rate in terms an integral 
over the product of the imaginary part of the electronic, 
dynamic susceptibility, and a form factor that encapsu- 
lates the interaction between the electronic and nuclear 



J^a 






_JSha 






J^ca 




A"" 



4 = 

and a by TT rotation around the c-axis, 
A. = 



j^aa J\^^ J^^*^ 

^ba jX^^ 



(19) 



(20) 



These four tensors can be combined using Eq. ([TT]) to 
find, 



A 



where, 



Aaa ^ ^ \ab^ ^ ^ Aac ^ ^ 

^ ^a^b ^ ^a^b ^a^b 

4 I -A'"'SaSi A'^^'CaCb iA'"'CaSb 

SaCb iA CdS}) A CaC}) 



(21) 



QadO 



C5 = COS 

Sb = sin 



qbbp 
2 

Qbbo 
2 



(22) 
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The internal field follows from substituting the above 
expressions for the nuclear-electron coupling tensors into 
Eq. dl]). We find, in agreement with 
internal field along the c-direction given by, 



12|, an average 





(hint) = ±4^"^(m^) I 



(23) 



where (m°) is the component of the average, electronic 
magnetic moment ahgned with the a-axis. The sign of the 
field depends on which As nucleus is under investigation. 




^ boqb 



aoQa 



FIG. 4: (Color online). q-dependence of the form factor 
J^(q, /igxt^), given in Eq. (|24|) . for an external field applied in the 
[110] -direction. The form factor acts as a "filter" of the electronic 
fiuctuations. It is finite at the ordering vector q = (Tr/ao, 0, tt/cq), 
and therefore allows the dominant fiuctuations of the longitudi- 
nal susceptibility to "pass the filter". We use the parameters 
^aa ^ Q QQ T/iiB and = 0.43 T / ijlb from O and make the 
approximation A^^ ~ A!^^ . 



Eq. (p!5]) and Eq. (|2T]) can be used to calculate the form 
factor for field applied in the [110]-direction as. 



[TlpiB) 




aoQa 



FIG. 5: (Color online), q-dependence of the form factor J-'(q, /^ext)? 
given in Eq. (|25|) . for an external field applied in the c-direction. 
The form factor is zero at q = (7r/ao, 0, gc), and therefore it "filters 
out" the dominant electronic fiuctuations at q = (7r/ao, 0, tt/cq). It 
has a maximum at q = (0, 0, 0), which matches the secondary peak 
in the imaginary part of the longitudinal susceptibility. We use the 
parameter A"'^ = 0.66 T/i^b from [IJ] and make the approxima- 
tion ^ 



electron coupling tensor is given by. 
1 



Mo 
47rrp 



-3af/r2 
-Zciai/rj 



1 



-36,Ci/rf 
l-3cf/rf 



(26) 



where = (a^, 6^, q) is a vector connecting the zth elec- 
tron moment to the nuclear site. This is longer range 
than the hyperfine coupling, but the symmetry of the 
nuclear environment remains C2v In consequence the 
position of the peaks and troughs in the form factor are 
unchanged, and thus the qualitative structure of the re- 
laxation rate will be the same. Since this form of cou- 
pling has been shown to be negligible in BaFe2As2 [ill, 
we concentrate exclusively on the hyperfine interaction. 



(24) 

where, following Eq. (j22j), c^, C5, Sa and 55 are q- 
dependent. For this field orientation, the form factor 
is finite at the ordering vector q = Q = (7r/ao, 0, tt/cq), 
as shown in Fig. (j4|). Exactly at this point 
J^(Q, /igxt^) = 16(^^^)^, and this is the approximation 
used in reference [19]. 

For field applied in the c-direction Eq. (p!6|) and 
Eq. (f2T]) can be combined to find, 

^(q, hl^,) = IQ{A^-Cac^f + lQ{A'^SaS^)\ (25) 

This has a peak at q = (0, 0, 0) and is zero at 
q= (7r/ao, 0, tt/co). Thus fiuctuations of the electron sys- 
tem at the ordering vector Q will be "filtered out" by the 
form factor (see Fig. (j5j)). 

In principle one could also include a dipolar coupling to 
the surrounding electron moments. The relevant nuclear- 



B. Dynamical, longitudinal susceptibility 

We now turn to the dynamical susceptibility of the 
electron moments, concentrating on the longitudinal fiuc- 
tuations relevant to BaFe2As2. The low energy field the- 
ory that reproduces the dispersion relation and has all 
the correct symmetries of the ordered state i o-'^^i^^i^^ — , 

^[^] = I d'rdr[h\^{drnf + ^ Pa{d^nf 

a—a,b,c 

-X±A^nl], (27) 

where x± is the static perpendicular susceptibility and pa 
is the spin stiffness along the ath crystallographic direc- 
tion. The relation between spin stiffness and spin wave 
velocity is Va = \/ pa/x^- The action is based on the 
non-linear sigma model, whose non-linearity arises from 
the requirement that = 1 in the partition function, 

Z = [ Vn6{n^ -1) e'^^^'l (28) 
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The constraint, n = 1, remains valid in the anisotropic 
model, Eq. (|27|) . and the anisotropic term, x±^^^a' 
arises from the fact that the a-axis is energetically 
favoured as the ordering direction in BaFe2As2. While 
the correct microscopic model for electronic magnetism 
in the pnictide materials remains controversial^—, we 
stress that this field theory provides a correct descrip- 
tion of their low-energy spin- wave excitations, regardless 
of the details of the high-energy physics. 

The longitudinal, dynamic susceptibility follows from 
the action given in Eq. (|27|) . In Appendix A we de- 
rive an expression for the susceptibility, using a Gaus- 
sian approximation to describe fluctuations of the order- 
parameter field, n, around the ordered state. This has 
two main contributions, one from q ^ Q and the other 
from q 0, and can be expressed as, 

^^{x||(q,^o)} ^^rn{x\\,sM ^ Q^^o)} 

+ {X||,nn(q ^ 0,cJo)} . (29) 

Taking the hmit cjq ^ in Eq. (IXT5]) and Eq. (|XT9]) 

gives, 

^^{X||,8t(q - Q^^o)} - ^^Y^^^^_^o 1 



(2^)^ 



'^k nB(a;i,k)(nB(^i,k) + 1) 



kBTxi 



L 



d^k nB(w2,k)(^^B(w2,k) + 1) 



(fi^2,k)2 



(30) 



and, 



'{X||,nn(q - 0,CJo)} ^ {gi/J^BfTrhVcell 



k.o (2^)^ 
(fk 



nB(^i,k)(nB(^i,k) + l)(^(^a;i,k - ^cji,k+q) 



of these peaks is controlled by the temperature, T, and, 
for all other wavevectors, the susceptibility is exponen- 
tially suppressed. For a realistic set of parameters the 
peak at q = Q is three of orders of magnitude larger 
than that at q = 0, as illustrated in Fig. ([6]). 



^^rn{x(q)} 

(10-1^ X ^l/s)^ 





h^qh 



(a) 



tT»"»{x(q)} 



(10-1^ X mIA) 




6oQ6 



f d k \ 

+ / 77r-^'nB{0J2M){f^B{0J2M) + l)(5(fiW2,k+q - fii^2,k) 

JkRiQ [^T^r ) 

(31) 

where S is the average electron spin per Fe-site, gi is the 
lande g-factor, is the standard Bose factor and. 



FIG. 6: (Color online). The q dependence of the imaginary 
part of the longitudinal, dynamic susceptibility at a) qc = tt/cq 
and b) qc = 0, as predicted by Eq. (|30|) and Eq. ((3T]) . The 
peak at q = (7r/ao,0,gc) is approximately 1000 times larger than 
that at q = 0. We use the parameters Va = vij = 280 meVA, 
Vc = 57 meVA and A = 9.8 meV from O and x± = 1-2 X 10-^ 
from The temperature is ksT = 1 meV. 



C. The relaxation rate with field in the 
[110]-direction 

We are now in a position to determine the relaxation 
' rate for an external field applied in the [110]-direction. 
The form factor is non-zero at q = Q = (yr/ao, 0, tt/cq), 
which is the peak in the imaginary part of the susceptibil- 
ity. Thus the neighbourhood of this point in momentum 
space will dominate the integral for the relaxation rate. 

Expanding the form factor in Eq. (|24|) around q = Q 
leads to. 



^(q 



Q,0 



16(^"°)2. 



(33) 



Sa;2,q= ./A2 + ^«2(q„-Q„)2 



(32) 



Substituting this into Eq. ([Tl]), along with Eq. (|3Q|) . and 
making the coordinate transformation. 



The imaginary part of the susceptibility has a large 
peak at the ordering vector q = Q = (7r/ao, 0, tt/cq) and 
a smaller peak at q = 0, as shown in Fig. (|6|). The width 



kla ^a^a 



(34) 
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gives, 

f (fkl f (fk2 1 / N/ / N N 

" Le Le (^"^('^'^^)(^^("'^^) + 

X S{huu\^^ — hcAj\^^)^ (35) 



where mo = QiIJ^bS, cjk = ^i,k, 



(36) 



is the geometric mean of the spin wave velocities and the 
integrals are over a cone of spin wave excitations. 

The density of states, which is given in Appendix [Bl 
can be used to transform the integral over momentum 
into one over energy, resulting in. 



1 hmlj%{aoboCo) 



ca\2 



X / de(e - A"') tt. 



(37) 



The required energy integrals are evaluated in Ap- 
pendix [CI and it follows that the relaxation rate is. 



1 f /[ca\2 f 



where. 



and. 



2fimg7^(ao6oCo )^A^ 



CstA 



(38) 



(39) 



^stAx) = x^Lhie-'^-) + x^Li2(e-i/^), (40) 
with Lim{z) = Xl^o ^V^^ polylogarithm of z. 

D. The relaxation rate with field in the c-direction 

We now turn to the relaxation rate with the external 
field applied in the c-direction. The form factor is qual- 
itatively different from that with the field applied in the 
[110]-direction, since it is no longer peaked at the order- 
ing vector but is in fact zero at this point. Expanding 
the form factor to lowest order around the wavevectors 
q = Q and q = gives, 

^(q « Q, « AiAn^aliQa - Qaf 

^(q«0,/i,^J«16(^-)2. (41) 

There are thus two main contributions to the relax- 
ation rate. The first, from the region around q = 0, we 



denote as 1/T{^^ and use Eq. (jSTJ) to write it as. 



J cone J CO 



1 



(27r)3 y,„„, (27r)3 {hu^^^f 

X (^(flCJki — ^^k2)- 

where the coordinate transformation. 



nB(c^ki)(nB(c^ki) + 1) 



(42) 



(43) 



has been applied. Making use of the results in Ap- 
pendix [B] and Appendix [Cl leads to. 



(A^^) Cun,l ^nn,l 



where, 



and, 



Co, 



8{gnB?h^%{aoboCofA'> 



(44) 



(45) 



^unA^) = a;2Lii(e-i/^) + 5x3Li2(e-i/^) 

+ 12a;^Li3(e-i/^) + 12a;^Li4(e-i/^). (46) 

The second contribution to the relaxation rate is from 
the region surrounding q = Q. This is suppressed rel- 
ative to the [110] field direction by the vanishing of the 
form factor at this point. We denote this contribution as 
1/Tf and find, 



Trhml^%{aoboCof 



1 

2 [' ■ \v, 
nB(wki)(nB(wki) + 1) 



(2^ 



cone 
2 



(27r)3 



26 + ^16 ) 



(47) 



where the coordinate transformation defined in Eq. (|34|) 
has been applied. Making use of the spectral functions 
derived in Appendix [B] and the integrals evaluated in 
Appendix O leads to a relaxation rate. 



where. 



st,k^ 



ba\2-\ 



^3d 



4hml-fj^{aoboCo )^A^ 



C, 



st,k'^ 



(48) 



(49) 



and. 



+ 3a;5Li4(e-i/^). 



(50) 
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The total relaxation rate for external field in the c- 
direction is given by the sum of the two contributions, 



Ti{h%,,) Tr{h%.,) Ti\K,,y ^ ' 

E. Comparison with experiment 

Having derived theoretical predictions for the relax- 
ation rate, we compare to experimental data, and show 
that the two are consistent at a quantitative level. 

The experimental data for BaFe2As2 are separable into 
an isotropic term, which is linear in temperature, and an 
activated anisotropic term, which has a more complicated 
temperature dependence. The isotropic contribution to 
the relaxation rate is likely due to a fluid of conduction 
electrons associated with ungapped portions of the Fermi 
surface. The linear temperature dependence would then 
be attributable to a Korringa-type relaxation rate^'^. 

We have previously argued that, at low tempera- 
tures, the interaction between these two electron fluids 
is negligiblei^. However, the origin of this isotropic re- 
laxation rate is not pertinent to this paper. We consider 
that the anisotropic term comes from the scattering of 
thermally excited spin waves and model the relaxation 
rate using the results of Sections IHI CI and IHIDi 

Describing the low temperature region of the exper- 
imental data with a function 1/Ti ^ CincT^ results in 
good fits to both data sets with Cine = 0.0325~^i^~^. In 
Fig. ([7j) this isotropic term is subtracted from the data 
and we concentrate on fitting the anisotropic contribution 
to the relaxation rate using the theory described above. 




FIG. 7: (Color online). Fits to the NMR relaxation rate data 
shown in the inset of Fig. (|2]). For external field in the [110]- 
direction (red circles) we fit Eq. p8|) . while for external field parallel 
to the c-axis (blue squares) we fit Eq. (I51|) . The same, linear T, 
isotropic term has been fitted and subtracted from both data sets. 
In fitting the anisotropic contribution to 1 /Ti we allow a single free 
parameter for each external field orientation: for the [110]-direction 
(A^^)'^Cst,i; and for the c-direction (A"'"')'^Cun,i- This provides 
convincing fits to the data, and we show in the text that estimating 
these fit parameters from independent experiments leads to quan- 
titative agreement. We use A = 9.8 meV, taken from neutron 
scattering experiments-^. 



Quantity 


Value 


ivietnoa 


iiei 


A 


y.o(^4j mev 


Neutrons 


1 ioj 


x± 


iu emu/mol — i.z x iu mev 


ivnignt sniit 


ri «i 
1^ 


[ao, bo, Co] 


To QH O TQ AT] /\ 
[Z.oU, Z. < y, D.4 < J /± 


Neutrons 




mo 


yj.oi fiB 


Neutrons 


ri /ii 


Aca 


n /I Q T' / 1 1 „ 


IN IVlxx 


ri 91 


Aaa 
J\ 


n fKP. T' / 1 1 ^ 

U.OD i / flB 


Knight shift 


ri 91 


7/v 


4t. D X lU ± S 


±aDuiatea 




Va 


280(150) meVA 


Neutrons 


[16] 


Vb 


280(150) meVA 


Neutrons 


[16] 


Vc 


57(7) meVA 


Neutrons 




hext 


1.5T 


NMR 


[12] 



TABLE I: Parameters used to fit the 1/Ti relaxation rate in 
BaFe2As2. 

For external field in the [110]-direction we treat 
(^^^)^Cst,i as a free parameter. This results in a con- 
vincing fit to the data for, 

{A'^fCst,i^7-^s-\ (52) 

Eq. (|39|) gives Cgt^i in terms of parameters that have been 
measured in independent experiments. Substituting in 
the values in Table |I] implies, 

0.13 < < 31 (53) 

where the error is predominantly due to the uncertainty 
in the spin- wave velocities. Therefore we find that the 
fits to the NMR data are in quantitative agreement with 
independent experiments, within the limits set by exper- 
imental error on measurements of the input parameters 
of the theory. 

We have previously shown that Eq. (j38j) also provides 
convincing fits to NMR 1/Ti relaxation rate data for 
SrFe2As2, with magnetic field parallel to [110] [l9|. Sim- 
ilarly, Klanjsek et al. have found good agreement with 
1/Ti data for NaFeAs 

For external field in the c-direction we fix the ratio, 

and in the absence of other information assume 
jaa ^ j^ha^ showu in Fig. ©, fitting Eq. ^ to the 
data with the single free parameter {A^^)'^Cun,i gives a 
convincing fit to the data for, 

0.12 < {A'^'fCun^i < 0.5 s-^, (55) 

where the uncertainty comes from the spin wave- 
velocities Va and V}). Estimating Cun,i using the values 
in Table H] gives, 

0.002 < (^^^)^Cnn,i < 0.25 s-^. (56) 

Again, within the bounds of experimental error, we find 
quantitative agreement between theory and experiment. 
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In summary, the theory for 1/Ti developed through- 
out this section is in quantitative agreement with the ex- 
perimental data of Kitagawa et a!^'^. This demonstrates 
the importance of taking angular resolution into account 
when calculating relaxation rates. We will now go on 
to generalise these results for an arbitrary strength and 
direction of the external magnetic field. 



IV. EXTENSION OF THEORY TO ARBITRARY 
ORIENTATION AND MAGNITUDE OF 
APPLIED FIELD 

In this section we develop a theory of the relaxation 
rate for arbitrary orientation and magnitude of the ex- 
ternal magnetic field. This follows from determining the 
general expression for the form factor in BaFe2As2, and 
then combining this with the above calculation of the 
longitudinal susceptibility. 

In order to find the form factor for an arbitrary mag- 
nitude and orientation of magnetic field, it is necessary 
to study the rotation matrices that transform between 
the (x, z) coordinates [those in which z is aligned with 
hnuc] and the (a, 6, c) coordinates [those aligned with the 
crystal axes]. Consider a rotation matrix R that ro- 

hext 

tates a vector from the (a, 6, c) coordinate system into the 
(x, z) coordinate system. The action of this matrix on 
the objects of interest is, 

a, (3 



a, 13 



(57) 



where a, /3 = {a, 6, c} and ^, = {x, z}. 

These rotation matrices can be used to transform 
Eq. (p!3|) for the relaxation rate into. 



1 



RV Rh 

next "ext 



Rr Rh 

next next 



.JO. .sp am{x"^(q,t^o)} 



hujQ 



(58) 



where 7, ^ = {a, 6, c}. 

Since only the longitudinal susceptibility, 
X^^(q, cjo) = X|| (q? ^o), is relevant to the relaxation 
process in BaFe2As2, it follows that, 

^^{xii(q. ^o)} 



_J - hm ^^y^.T V h r^^^ll^"^ 

^]Jl^) - .0^0 2N^^^ ^ ^^"^^ ^^^^^ h^o 

^^^^ «^T3A/rR7 ^ 



qGPMBZ 



(59) 



where the form factor that couples to the longitudinal 
spin fluctuations is. 



J'(q,hext) = [• 

7,(5 



RZ 



TDX6 

next 



RV 

nex1 



R 



y5 



6a 



(60) 



For an external field of arbitrary magnitude, the as- 
sumption hnuc ^ hext is uo longer valid. As such it is 
natural to define two sets of angles, as shown in Fig. (j8j). 
The first set of angles, ('0, A), describe the orientation of 
hnuc in the (a, 6, c) crystallographic coordinate system, i/j 
is the angle between the z-axis and the c-axis and A is 
the angle between the projection of the 2:-axis onto the 
ab plane and the a-axis. The rotation matrix is thus. 



^sin^ A + cos t/j cos^ A 
— sin2Asin^ ^ 
— cos A sin i/j 



-sin2Asin^ | 
cos^ A + cos t/j sin^ A sin A sin 



sin A sin ip 



cos A sin tjj 
1 A sin 



(61) 



The second set of angles, {9,4'), describe the orientation 
of 

hext ill the (a, 6, c) crystallographic coordinate system. 
These are the experimentally accessible set of angles, 
where is the angle between the c-axis and hext, and 
(j) is the angle between the a-axis and the projection of 
hext onto the a6-plane. 



h J 




h 






^ext 














^.-^^^ I'l' ^^^^^^^ 





FIG. 8: (Color online). The relationship between the coordinate 
system of the crystal axes of BaFe2As2, (a, 6, c), and the coordinate 
system of the effective magnetic field hnuc at the '^^As nucleus, 
(x, y, z). In the magnetically ordered phase, the electron moments, 
m, are orientated along the crystallographic a-axis. The interaction 
between the ^^As nucleus and these magnetically ordered electrons 
gives rise to an effective internal magnetic field, hjnt, directed along 
the c-axis. The total field hnuc = hext + hjnt is the sum of hjnt and 
the external magnetic field, hext, applied during NMR experiments. 
The orientation of the external field, hext, relative to the crystal 
axes [shown here with polar angles (^, (j))] can be varied at will 
by rotating the sample in a goniometer. This in turn changes the 
orientation of the total effective field hnuc [shown here with polar 
angles ('0, A)]. 

For the As nucleus in BaFe2As2, Eq. ([23|) gives. 
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(hint) = (0,0,±|(hint)|), and it follows that, 



to find, 



x = <p, 

tan?/^ — 



I hext I sin l9 



hext|cOs6>± I (hint)!' 



(62) 



Thus the {ip, A) angles that enter the theory can be ex- 
pressed in terms of the known angles {0,(t>). In the 
high external field regime, |hext| ^ |(hint)|, and there- 
fore ip = and A = 0. 

The form factor that follows from substituting Eq. ()UT|) 
into Eq. (|60| is, 



7"(q, hext) = (cos^ + sin2 A sin2 ^p) A^^A"l^^ 

+ (cos^ + cos^ A sin^ A'^^A''^^ 
+ sin^ V A^^A"^^ 

- i sin 2A sin^ ^ (^^M^L^q + A'^^^A 
+ ^ cos A sin 2i> + 



sin A sin 2tP (^^Mi"q + A^^A^) . 



(63) 



Eq. (|2T]) can be used to re-express this as. 



J'(q, hext) = 16 (cos^ + sin^ A sin^ V^) (^^^c^c^)^ 
+ 16 (cos^ + cos^ A sin^ ^) {A^^'saSbf 

+ 16sin2Asin^ ^^M^^CaC55aS5, (64) 

where c^, c^, 5^ and are given in Eq. (|22]) . Approxi- 
mating the form factor close to q = Q gives, 

J'(q ^ Q,hext) ^ 

16(^"^)^sin2V^ 

+ 4(^^^)2 (cos^ ^ + sin^ Asin^ V^) ag(ga - Qaf 
+ 4(^^'^)2 (cos^ i/j + cos^ Asin^ V^) ^^(g^ - Qbf 

+ 4^^M^"sin2Asin2 ao6o(ga - Qa)(g6 - Qb), 

(65) 

while close to q = the leading contribution is. 



^•(q ^ 0, hext) ^ 16(^^^)^ (cos^ t/j + sin^ A sin^ V^) . 



(66) 



The general form of the relaxation rate is now accessi- 
ble. The techniques outlined in Section [Till can be used 



+ {A"")^ (cos^ V + sin^ Asin^ V) C„„,i <E>„„,i 



(cos^ ^/^ + sin^ A sin^ t/j) 



+ (cos2 ^ + cos2 Asin^ ^) ^^^^] ^st,k^ . 

(67) 



This equation for 1/Ti leads to a "doughnut" -shaped an- 
gular dependence, as illustrated in Fig. (|9]). The rate is 
largest when hext is orientated in the a6-plane, and small- 
est for hext in the c-direction. There is a small difference 
between external field aligned in the a-direction and in 
the 6-direction, with the 6-direction being faster. 




'*ext 



FIG. 9: (Color online). "Doughnut" shaped prediction for the 
variation of the ^^As NMR 1/Ti relaxation rate of BaFe2As2 with 
the orientation of the external magnetic field hext, at fixed temper- 
ature. The value of 1/Ti is represented by the radial distance of 
the surface from the origin, and calculated from Eq. ()67|) using the 
parameters given in Table [l] for T = 0.6A, |hext| = 6T. The relax- 
ation rate is a minimum when the external field is applied along 
the c-axis, parallel to the internal field (hint). 



As shown in Fig. ^TO\i , the relaxation rate depends not 
only on the orientation of the external magnetic field, but 
also on the magnitude. This arises from the fact that it 
is the orientation of hpuc that determines the relaxation 
rate. When hext = 0, it follows that hpuc = hjpt, and the 
nuclear field is parallel to the c-axis. If the external field 
is then applied along the C-axis, the orientation of hpy^ 
remains unchanged. If hext is slowly turned on in the 
a6-plane then, as the magnitude is increased, hnuc will 
rotate by 7r/2 radians. Thus different strengths of the 
internal field will correspond to different angles (?/^, A) of 
hnuc, even when the (^, </>) angles of the external field are 
kept constant. 

For an external field applied in the 6-direction the re- 
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FIG. 10: (Color online). Theoretical predictions for the 
anisotropic contribution to the relaxation rate as the magnitude of 
the external field, hgxt, is varied. For an external field applied in the 
c-direction there is no dependence on magnitude. For an external 
field applied in the 6-direction the relaxation rate is strongly depen- 
dent on magnitude, as described by Eq. (|68|) . We set T = 0.6 A and 
use the parameters obtained from fits to BaFe2As2 data (Eq. (|52|) 
and Eq. (|55]) 1. 



laxation rate is given by, 

u2( Aha\2 



Mlr^r * ^^^^ 



A ; Khint)|2 + |hextP 



(^°°) Cun,l ^un,l 



knT 



2 ^st,k^ '^st,k^ 



(68) 



where we emphasise the dependence on the magnitude of 
hext, and the constants Cst,i, C„„4 and Cst,k^ are defined 
in Eq. Eq. (gSI and Eq. If the external field 

is applied in the a-direction then, 

1 ^,(Aca-.2f^ ^ ( kBT \ IhextP 



V A ; Khint)|2 + |hextP 



5 ^S^,fe2 ^t>s^,fe 



A J |(hint)P + |hextP 

2 ^St,/C2 ^St,/C2 — ^ . (69) 



In both cases the first term is the dominant contribution. 

An experimental measurement that either mapped out 
an octant of the "doughnut" shown in Fig. (|9]), or one 
that showed the dependence on external field magnitude 
illustrated in Fig. (fTQ|) . would provide strong support for 
the angle-resolved theory developed in this paper. 



V. QUANTITATIVE DETERMINATION OF 
SPIN WAVE VELOCITIES FROM NMR 

Having established that it is possible to quantitatively 
match the theory of the NMR relaxation rate to ex- 
periment, we show how this can be exploited to make 
quantitative measurements of the dynamical properties 
of collinear anitferromagnets, with BaFe2As2 as an ex- 
ample. This technique could prove especially useful when 
crystal sizes are too small for inelastic neutron scattering 
measurements, as is often the case for newly synthesised 
materials. 

NMR measurements of the internal field provide an 
excellent way of determining static properties of antifer- 
romagnets, such as the ordering vector, the direction of 
the easy-axis and the energy scale of the resultant gap 
in the spin-wave spectrum. We have shown previously^^ 
how the gap. A, can be measured in this way. Also, 
Knight shift measurements allow the components of the 
nuclear-electron coupling tensor to be determined^^. 

We now argue that NMR relaxation rate measurements 
can be used to determine the hydrodynamic properties 
of collinear magnets, and in particular their low-energy 
spin wave velocities. We note that it has previously been 
suggested that it might be possible to place bounds on 
spin wave velocities in magnetic Fe pnictides from NMR 
experiment^. However we believe that ours is the first 
theory sufficiently advanced to make a quantitative com- 
parison with experiment. 

The first step is to determine the geometric combina- 
tion = VaVhVc by fitting data for the relaxation rate 
with external field in the a6-plane using Eq. (|38]) . 

Next, we consider measuring the relaxation rate along 
the three principal crystal axes, and combining these in 
the linear combinations. 



1 



1 



± 



1 



1 



T± 2 [T^{hl,,) 1 - sin^ e V^l(^ext) Ti(/le^t) 

(70) 

The form factors associated with these combinations are 
constructed from Eq. ([64|) . and given by, 

^+(q) = lQ{A'^SaS,)\ 



in the case of and, 

J^-(q) = 16(^""CaC6)2, 



(71) 

(72) 



for , where c^, C5, Sa and are given in Eq. ([22]) . 

It follows that, 



C, 



&0 



V A y 



and, 



^st,k^ I 

1 . fkeT 



(73) 



1^ \ A 



^st,k- ( ^ ) , (74) 
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where |hext| ^ I (hint) I is assumed and the constants Cst,i^ 
Cun,i and Cst,k'^ are defined in Eq. (|39|) , Eq. (|^5]) and 
Eq. (|49|) . The only remaining unknowns are the spin 
wave velocities Va and and therefore these can be ex- 
tracted from experimental data via, 



1/T+ 



and, 



1/Tf - (^-)2c„„,i$„„,r 



(75) 



(76) 



The velocity Vc then follows from the value of ^5. 

In order to perform a check on the values of the spin 
wave velocities one could also measure the dependence of 
l/Ti(/igxt) and 1/Ti{hl^^) on |hext| at constant tempera- 
ture, as illustrated in Fig. Eq. (|68j) and Eq. §9^ 
can be combined to give. 



1 



Cst,k^^st,k^ 



bl{A^ [ksT 
vl 

IhintP 



|(hi„t)|2 + |hextP^ 



(77) 



where Va and Vh are the only free parameters. 

These techniques for measuring individual spin wave 
velocities involve combining relaxation rate measure- 
ments such that the dominant processes are canceled out, 
and the subleading terms are revealed. Thus a high de- 
gree of experimental accuracy is required. However, if a 
nuclear site can be found at which the internal field van- 
ishes, then the fiuctuations of the electron moments at 
the ordering vector are "filtered out" by the form factor 
for all orientations of internal field. In this case the above 
techniques are likely to become more powerful, since the 
cancellations required will be smaller. For example, this 
appears to be the case for the Y nucleus in YBa2Cu306 



VI. CONCLUSION 

In this paper we presented a theory of NMR 1/Ti re- 
laxation rates in a collinear antiferromagnet that pro- 
vides quantitative fits to published data for ''^As NMR 
in BaFe2As2. All predictions are given in absolute units, 
and the spin fiuctuations of electrons are parameterized 
in terms of a small number of hydrodynamic parameters 
— the ordered moment mo, transverse susceptibility 
anisotropy gap A and spin wave velocities {va^vi)^Vc). 
The remaining parameters of the theory are the small 



number of matrix elements of the transferred hyperfine 
interaction between Fe electrons and the ^^As nuclear 
spin. Since these can be determined from measurements 
of NMR spectra, the resulting theory has no adjustable 
parameters. 

A key feature of this theory — and of ^^As NMR on 
BaFe2As2 — is a strong dependence of the 1/Ti relax- 
ation rate on the orientation of the magnetic field. This 
angle-dependence can be traced back to the "filtering" of 
spin fiuctuations by the form factor for transferred hy- 
perfine interactions, which in turn depends on the orien- 
tation of the magnetic field. Taking this into account, the 
theory correctly captures the qualitatively different tem- 
perature dependences of 1/Ti for ^^As NMR in BaFe2As2 
with field applied along the [110] and [001] direction^. 

Moreover, since the theory is expressed only in terms 
of hydrodynamic parameters of the magnetic electrons, 
this fitting procedure can be inverted, and angle-resolved 
1/Ti measurements used to determine spin wave veloc- 
ities directly from NMR experiments. We have pro- 
posed a specific scheme for doing this from ^^As NMR 
in BaFe2As2. 

While we have developed this theory with the par- 
ticular goal of explaining ^^As NMR experiments in 
BaFe2As2, the results have a much wider applicability. 
Firstly, a similar analysis can be applied to other collinear 
magnets, simply by modifying the form factor to take 
into account the symmetry environment of the nucleus in 
question. It appears that the "doughnut" shaped angle- 
dependence of 1/Ti, shown in Fig. (j9j), remains valid for 
all nuclei that experience a non-zero internal magnetic 
field hint- The "hole" of the "doughnut" is aligned with 
hint- For nuclei at high- symmetry sites where this inter- 
nal field vanishes, the leading term in the relaxation rate 
is "filtered out" for all directions of external field, and the 
angular resolution acquires a more isotropic "peanut" - 
like shape^^-. 

Although the theory developed in this paper is specific 
to an ordered antiferromagnet, the idea of angular resolu- 
tion in 1/Ti measurements can easily be extended to the 
study of critical fiuctuations. Relaxation rate data for 
BaFe2As2 with field applied in the [110]-direction show 
a significant upturn in 1/Ti as the (first-order) magnetic 
phase transition at Tn = 135 i^^ is approached from the 
paramagnet [cf. Fig. ([1])]. This upturn occurs because, 
for field parallel to [110], 1/Ti probes spin fiuctuations 
near to the magnetic ordering vector q = Q, and these 
are enhanced approaching the phase transition. In con- 
trast, when the field is applied in the [001] direction, 
there is no upturn in 1/Ti. For this field orientation, the 
form factor "filters out" critical fiuctuations at q ^ Q, 
and 1/Ti is determined instead by spin fiuctuations with 
q ^ 0. Angle-resolved NMR experiments can therefore 
be used to isolate critical fiuctuations near to a phase 
transition in BaFe2As2 and other antiferromagnet^^. 

It can be seen in Fig. ([1]) that the relaxation rate re- 
tains a significant angular dependence in the paramag- 
netic phase. For an appreciable range of temperatures 
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above T/v the susceptibility will be dominated by fluc- 
tuations close to the ordering vector. Thus the idea of 
angle dependent "flltering" of spin fluctuations by the 
form factor will remain important. This is a theme that 
we will develop in a separate publication^^. If the tem- 
perature is further increased such that T ^ Tn^A then 
the susceptibility will become isotropic in spin space and 
only weakly q-dependent. In this regime the angular res- 
olution of the relaxation rate will be primarily due to the 
anisotropy of the nuclear-electron coupling tensor^^. 

This ability to tune between different spin fluctuations 
should also make angle-resolved NMR a powerful probe 
of unconventional magnetism, and in particular of ex- 
otic quantum phases in frustrated magnets. The absence 
of an ordered magnetic moment, lack of a large single 
crystal, or the requirement of large magnetic flelds, often 
make these systems inaccessible to other probes, such 
as neutron scattering. One intriguing possibility is that 
angle-resolved 1/Ti measurements could provide a posi- 
tive means of identifying the long-sought quantum spin- 
nematic state, a magnetic analogue of a liquid crystal, 
which does not break time reversal symmetry, and so 
does not give rise to magnetic Bragg peaks or static split- 
ting in NMR spectrg^ — . This theme will be developed 
elsewhere^. 
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Appendix A: Calculation of the longitudinal 
suceptibility 

In this Appendix we calculate, at Gaussian order, the 
imaginary part of the dynamic susceptibility due to lon- 
gitudinal fluctuations of the ordered moments. While the 
calculation is not original^^, we are not aware of any pub- 
lished derivation in absolute units, and therefore include 
it here for completeness. To make the "flltering" effect of 
the form factor more transparent, we choose to work in 
the Brillouin Zone associated with the orthorhombic lat- 
tice of magnetic sites in the paramagnet (PMBZ), rather 
than the magnetic Brillouin Zone of the low temperature 
antiferromagnet . 

The dynamical susceptibility of the non-linear sigma 
model Eq. (|27|) can most easily be calculated by consid- 
ering the effect of an external fleld which varies in space 
and time, h(r, r) [meV]. The real space susceptibility 
[/i^] in the (a, 6, c) coordinate system of the crystal lat- 
tice is then given by. 



Sh"{r,T)5hl3{0) 



■ (Al) 



where F = —\nZ^ and Z is the partition function, which, 
for the non-linear sigma model, is given by, 

Z = Jvn S{n^ - 1) e-^["'^l (A2) 

The longitudinal susceptibility can be accessed by con- 
sidering both staggered and uniform flelds applied paral- 
lel to the ordering axis. We flrst consider the effect of a 
staggered fleld, hgt- This couples directly to the antifer- 
romagnetic order parameter, n, according to 



5[n,h5t] = S[n] 



hVr 



cell 



d rdr n.ho 



(A3) 



where S is the total spin per site. The order parameter 
fleld can be parametrised as n = 



1 



and, for temperatures at which the fluctuations around 
the ordered state are small, ^1,^2 <^ 1- To Gaussian 
order in (/)i,(/)2, Eq. (|27|) becomes. 



1 



2hVeeU 



(frdr 



Oi 

+X±AV'] , (A4) 



where = (0i,02)- Likewise Eq. (|A3[) for the staggered 
fleld gives. 



S[(P,hst]=S[4>\ 



Shs 



hVcell 

Starting from the partition function, 



(frdT ( 1 - — 



(A5) 



(A6) 



and calculating the longitudinal, staggered susceptibility 
from Eq. (|A1|) gives. 



X||,,(r,r)=ry(r) (^^^ ((^^(r, r)0^(O)) 

-{4>\v,r)){4>\0))), 

(A7) 

where 77 (r) accounts for the staggered nature of the fleld 
and. 



(O) = ^ j V(j)0 e-^[^l 



h=0 



(A8) 



The susceptibility in Eq. p3|) is the Fourier transform 
of that entering Eq. (|A7[) . We choose to work in the full 
PMBZ in reciprocal space, rather than the reduced mag- 
netic Brillouin Zone (MBZ), since this makes the physi- 
cal picture of the interaction with the form factor clearer. 
For two sublattice antiferromagnetic order there are two 
identical cones of spin wave excitations in the PMBZ, 
one at q = and the other at the ordering vector q = Q. 
After Fourier transform, the flelds ^1 and 02 describe 
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two independent cones of bosonic excitation. The field 
Fourier transforms are defined by, 

iuJri 
iuJri 

(A9) 

where (j)i{q_^iujn) describes the excitation cone at q = 
and 02(q, ^^n) describes the excitation cone at q = Q. 
It follows from substituting the field Fourier transforms 
into the action that the two field averages are, 

(0i(q,ia;n)^i(q ,^^n)) = 

^(q + q!)S{iuJn + iuJn)Gp (q, iuJn) 

(^2(q,^CJn)(/>2(q',^^n)) = 

^^(q + q - ^Q)S{iuJn + ioj'JGp (q, iuJn) 



6i(q,ia;n)^2(q ,^^n)) = 0. 



(AlO) 



with. 



Gp{q,iuJn) = 2 -r- ^ 



The energies. 



(All) 



(A12) 




describe the dispersion of the two spin wave cones and. 



{A13) 



gives the spin wave velocities [meVA], where 
a = {a, 6, c}. 

Expanding four field averages using Wick's theorem, 
using Eq. (jAlOp to substitute for the two field averages, 
Fourier transforming the fields using Eq. (|A9p and rewrit- 
ing the staggering parameter as 77(r) = e'^*^ '^, gives, 

/ . X f gii^Bsy ^ 

X||,sHq,^^n) - ( ) Vcell 



\ 2 J 
1 



(A14) 



where gi is the Lande g-factor. Performing the Matsub- 
ara sums over iun and analytically continuing to real fre- 
quencies we find 



^^{X||,5t(q,^0)} ^ (^^^^^^^ ^^"^cell-^ 

( f d^k nB(c^i,k) - nB(c^2,k+q) 

X [(5(flCJl,k - ^CJ2,k+q + ^O) - (^(^^2,k+q " ^^l,k + ^o) 

d^/c nB(a;2,k) - ^B(^i,k+q) 



k~Q (27r)3 fia;2,k ^cji,k+q 

[(5(flCJ2,k - ^^l,k+q + ^O) - (^(^^l,k+q - /^CJ2,k + ^o)] ] , 

(A15) 



where cji^q = C(j2,q+Q has been used and. 



^^(^q) 



(A16) 



is the standard Bose factor. The staggered susceptibility 
describes scattering of spin waves between two cones of 
excitations separated by wavevector q = Q. It is peaked 
at q = Q, and the sharpness of the peak increases with 
decreasing temperature. This is because, at low tem- 
peratures, less of the excitation cone is accessible to the 
spin- wave fluctuations. 

There is also a contribution to the susceptibility asso- 
ciated with the application of a uniform field, hun- The 
couphng between the order parameter field and hun can 
be included in the action a^*^, 



S[n,hun] = S[n] - — ^ / d^rdr [ix±hun'{n x 9^11) 

Vcell J 

+ ^ (hL - (n.h„„)2) - — ^ada{n.hu„) 

(A17) 



Performing an expansion in the fields 0, as before, leads 
to 



S[(l),hun] = 



cell 



d'^rdr 



un fi2 I /2 



2h 



(A18) 



After repeating a similar set of manipulations to above, 
the contribution to the imaginary part of the susceptibil- 
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ity is found to be, 



'{X||,nn(q,^o)} ^ {gil-lBf-hVcell 



[L 



3 (nB(c^i,k) - nB(c^i,k+q)) 



X [5{hWi^\^ - fia;i,k+q + ^o) - (^(^^l,k+q - ^CJi,k + CJo)] 



kf^Q 



(27r)3 (^^(^2,k) - nB(cJ2,k+q)) 



X [(5(flCJ2,k - ^^2,k+q + CJo) - (^(^^2,k+q - ^^2,k + CJq)] ] 

(A19) 

This describes scattering of spin waves within an exci- 
tation cone, and is peaked at scattering vector q = 0. 
As with the staggered susceptibihty, the peak becomes 
sharper at lower temperatures, where less of the spin 
wave cone can be accessed. 

The imaginary part of the susceptibility is given by the 
sum of the two terms, 

^^{xii(q,^o)} =^^{xii,s^(q ^ Q,^o)} 

+ ^m{X||,nn(q^O,CJo)}. (A20) 

These results are used in Section UlIB I of the paper. 

1. Spin wave theory from Heisenberg model 

Since experiments on magnetic Fe pnictides are often 
discussed in terms of a Heisenberg mode l^^i^^ , we sketch 
below an equivalent calculation of the longitudinal sus- 
ceptibility within conventional spin wave theory. We 
stress that the non-linear sigma model correctly repro- 
duces the low energy behaviour of the Heisenberg antifer- 
romagnet, as it does for any microscopic model with the 
correct symmetries. We provide in Table |IT1 a dictionary 
to translate between these two models at the Gaussian 
level of approximation. 

The Heisenberg Hamiltonian with the correct symme- 
tries for the magnetically ordered phase of BaFe2As2 isi^, 

= Jia ^ Si.Sj + Ji5 ^ S^.Sj + Jic ^ S^.Sj 



(A21) 



where (ij) la counts first-neighbour bonds in the a- 
direction, {ij)2 second- neighbour bonds in the a-b plane, 
and Kab and Kc are single-ion anisotropics. 

We consider the case Kab = Kc = and use linear 
spin wave theory to calculate the spin wave dispersion 
and longitudinal susceptibility. Due to the large (A ^ 
100 K) gap in the spin- wave spectrum, 1/S corrections 
are small and can be safely neglected. By rewriting the 



spin degrees of freedom in terms of Holstein-Primakoff 
bosons^i one can transform the Hamiltonian to. 



Hisw = ^ (4'«-k) ( 

l^ePMBZ \ 



(A22) 



where. 



Ak = 2S{2J2 + Jia - Jib{l - coskb) + Jic)^^SK 

= 2S{2J2 cos ka cos kb + Jia COS kc + Jlc COS kc). 

(A23) 

Performing a Bogoliubov transformation with the coher- 
ence factors, 



^.k = -^J^ + 1, ^k = -^J^ - 1 (A24) 



V ^k 



V ^k 



results in. 



(A25) 



It follows that the imaginary part of the longitudinal sus- 
ceptibility is. 



^^{xii(q, ^o)} 



2 N 

^(iik^k+q-Q + ^k^k+q-q)^ (^^(^k) - ^^(c^k+q-Q)) 
kGPMBZ 

[S{huJ\^ - flCJk+q-Q + ^o) - ^(^^k+q-Q " ^^k + ^o)] • 

(A26) 

The combination of coherence factors appearing in the 
expression for the susceptibility can be rewritten in terms 
of the parameters of the non-linear sigma model accord- 
ing to, 

(^k^k+q-Q + ^k^k+q-q)^ ~ 

u;ku;k+q-Q (2I7) ' q ^ Q (A27) 
1, q^O 



with, 



(A28) 



The relationships between these hydrodynamic parame- 
ters and the exchange integrals of the Heisenberg Hamil- 
tonian are shown in Table |TT1 



Appendix B: Spectral representations of form factors 

The simplest way to perform the momentum integrals 
that occur in Section IIII CI and Section IIIIDI is to make 
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non-linear 


Heisenberg 


A 


4S[K(Jia + 2J2 + Jic) + K^]i 


Va 


2S[AJl - 4JiaJ2 + J?a + 2^2 Jlc - JlaJlc] ^ 


Vb 


2S[4J^ - 2JiaJ2 + 2Jib J2 + 2J2 Jlc 


Vc 


2S[2J2Jic- JiJic + JU^ 


X± 


l/(4Jia+8J2+4J^+4K) 



This is calculated as, 

dVt .0. 9 , /" „,4e(5(A;- Ve2-A2) 



(27r)- 



■ sin^ ^ cos^ 6 / d/c/c^ 



e(e2-A2)S. 



67r2 



By a completely analogous method, 



cone 

1 



(27r)= 



67r2 



e(e^-A^)t. 



Ve2 - A2 



h^{f^- ek) 



(B6) 



(B7) 



TABLE II: Relationship between the hydrodynamic parameters 
which characterize the non-linear sigma model and the exchange 
integrals which enter the Heisenberg Hamiltonian. Both models are 
treated at the same level of approximation : linear spin wave theory 
for the Heisenberg model and a Gaussian mean-field approximation 
for the non-linear sigma model. 



a transformation from momentum to energy space. This 
is done via the density of states, 



5(e) 



(2^ 



8{e - ek), 



with, 



ek = v/A2 + fc2. 



(Bl) 



(B2) 



Appendix C: Integrating products of Bose functions 

The calculation of the relaxation rate requires integra- 
tion of products of polynomials and Bose functions. We 
show here the results that we make use of in Section UlI CI 
and Section ITlIDi For integer n, 



/ 



da; 



(e^ - 1)2 



E 

m=0 



{n — m)\ 



where, 



(01) 



The integration region is spherically symmetric and so it 
is natural to use polar coordinates. Using, 



Ve2 - A2 



leads to a density of states. 



= ^ev^^. 



(B4) 



The expression for the relaxation rate with field par- 
allel to [001], Eq. (|17|), contains additional factors of /c^ 
and kl in the integrand. These are most easily handled 
using the spectral representation, 

^^"(^)=/ 7^^a'^(^-^k). (B5) 

J cone V^''/ 
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^=0 



The integrals required are. 
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